Ab initio derivation of electronic low-energy models for Ceo and aromatic compounds 
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We present a systematic study for understanding tiie relation between electronic correlation and 
superconductivity in Ceo and aromatic compounds. We derived, from first principles, extended 
Hubbard models for twelve compounds; fee KsCeo, RbaCeo, CsaCeo (with three different lattice 
constants), A15 CssCeo (with four different lattice constants), doped solid picene, coronene, and 
phenanthrene. We show that these compounds are strongly correlated and have a similar energy 
scale of the bandwidth and interaction parameters. However, they have a different trend in the 
relation between the strength of electronic correlation and superconducting transition temperature; 
while the Ceo compounds have a positive correlation, the aromatic compounds exhibit negative 
correlation. 

PACS numbers: 74.20.Pq, 74.70.Kn, 74.70.Wz 



I. INTRODUCTION 

Superconductivity in 7r-clectron systems, whose his- 
tory dates back to studies on graphite intercalation com- 
pounds in GO'Sfi has attracted broad interest in con- 
densed matter physics. Recently, two seminal discover- 
ies for carbon-based superconductors have been reported. 
One is A15/fcc CssCeo, with a new method to synthe- 
size highly crystalline samples<^"— The observed super- 
conducting transition temperatures (Tc) for the high- 
pressure samples are found to be as high as 38 K for A15 
and 35 K for fee. The other is potassium-doped solid 
picene;^ which opened a new avenue to various ' aromatic 
superconductors\ for which the maximum Tc has reached 
33 K^-ia 

The mechanism of superconductivity of these new su- 
perconductors has not been fully understood. For alkali- 
doped fuUeridcs, while there are several experimental 
reports which seemingly support the conventional BCS 
mechanism, various indications for unconventional super- 
conductivity have been also observed. For example, al- 
though the positive correlation between Tc and the lattice 
constant found in K- and Rb-doped fullerides has been 
understood in terms of the standard BCS theory^ more 
recent experiments for larger cations have revealed that 
Tc does not necessarily behave monotonically as a func- 
tion of the lattice parameters)^"— The fact that the super- 
conducting phase has a dome-like shape in the phase di- 
agram is indeed reminiscent of cupratesiSi and unconven- 
tional organic superconductors such as BEDT-TTF In 
addition, these new Ceo superconductors are insulators 
at ambient pressure?^"— indicating that the superconduct- 
ing phase resides in the vicinity of an insulating phase. 
In fact, considering these characteristic features, it has 
been proposed that interplay between orbital, spin, and 
lattice degrees of freedom is the origin of the high Tc 
superconductivity.— 

For aromatic superconductors, following the discovery 



of K-doped picene (C22lli4))^ it has been found that var- 
ious hydrocarbon compounds such as coronene (C24H12), 
phenanthrene (C14H10), and l,2:8,9-dibenzopentacene 
(C30H18) also exhibit superconductivity<i"— Diversity 
of hydrocarbon molecules suggests the possibility of 
new and higher Tc aromatic superconductors. So 
far, electronic structure)^"— electronic correlations ^^i^ 
electron-phonon interactions ^2*2^ and exciton/plasmon 
properties of the normal stateSir— have been studied, 
while studies on superconducting properties are still lim- 
ited. The pairing mechanism is totally an open question. 

There are several similarities between Cgo and aro- 
matic superconductors; they are both molecular solids 
having narrow bands around the Fermi level, whose en- 
ergy scale competes with that of electron-phonon and 
electron-electron interactions. The competition among 
these three factors is a characteristic aspect of carbon- 
based materials. Ah initio derivations of effective low- 
energy models for these compounds are important to 
make the situation transparent and to clarify the origin 
of their high Tc superconductivity. By comparing the 
parameters in the effective models for the Cgo and aro- 
matic superconductors, the differences and similarities 
are quantitatively identified and analyzed. 

Recent methodology for construction of the electronic 
model, based on the combination use of the maximally 
localized Wannier orbital (MLWO) (Ref. HI) and the 
constrained random phase approximation (cRPA);^ ex- 
tends its applicability range. It does not place limitations 
on the character of basis orbitals of the effective model, 
whether atomic or molecular. Indeed, it has already been 
applied to various complex systems such as BEDT-TTF 
(Refs. [13 and [31I ) or zeolites<^ While electronic inter- 
action parameters of the Cgo and aromatic supercon- 
ductors have been estimated by various methods )^^'^'^" — 
explicit and direct comparison of these systems by the 
same method has yet to be done. Thus it is imperative 
to evaluate the interaction parameters of Cgo and aro- 
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matic superconductors by exploiting the state-of-the-art 
technique and perforin a systematic comparison. In the 
present study, we constructed ab initio extended Hub- 
bard models which describe the low-energy electronic 
structure of twelve examples of Cgo and aromatic com- 
pounds. The transfer integrals were given as matrix el- 
ements of the Kohn-Sham Hamiltonian in the Wannier 
basis. The interaction parameters were evaluated by cal- 
culating the Wannier matrix elements of the screened 
Coulomb interaction, which is obtained by cRPA. The 
estimated correlation strength as the ratio of the interac- 
tion energy to the kinetic one is nearly or beyond unity 
for the studied materials, indicating that both the Ceo 
and aromatic systems are classified into a strongly corre- 
lated electron systems. On the other hand, we observed 
a notable difference between the two systems; for the 
Cgo system, there exist positive correlation regime in the 
correlation strength and the experimental Tc- In con- 
trast, the aromatic system exhibits negative correlation 
between these two quantities. 

This paper is organized as follows. In SeclU we show 
how to construct the low-energy models from ab initio 
calculations. In Sec lIIIl we show the calculated band 
structure, MLWOs, transfer integrals, and effective in- 
teraction parameters. We discuss the material depen- 
dence of the derived parameters and relation between the 
strength of electronic correlation and superconductivity 
in Sec lIVI Finally we give a summary in SecfVl 



II. METHODS 

We derive electronic low-energy models with the com- 
bination of MLWO and cRPA. This method has widely 
been applied to the derivation of effective models for 
3d transition metals,— their oxides,— organic conduc- 
tors^Si^ zeolites?^ iron-based superconductors )2M£ and 
5d transition metal oxides4i We first perform band cal- 
culations based on density functional theory (DFT) ,— 
and choose 'target bands' of the effective model. By 
constructing MLWOs for the target bands, we calculate 
transfer integrals and effective interactions in the effec- 
tive model. In the calculation of the effective interaction, 
the screening by electrons besides target-band electrons 
is considered within cRPA (see below). 

We apply this scheme to the derivation of effective 
models, i.e., extended multi-orbital Hubbard models, for 
the Ceo and aromatic compounds. The Hamiltonian con- 
sists of the transfer part Tit , the Coulomb repulsion part 
Tiu, and the exchange interactions and pair hopping part 
Hj defined as 

n = nt + nu + nj, (i) 



where 

Ht=^^^irim(Rii)a^^a^„, (2) 

(J ij nm 

t^-(R^.>^^<«;^<' (3) 

(Tp ij nm 
ap ij nrn 

with a^J (O'in) being a creation (annihilation) operator of 
an electron with spin a in the n-th MLWO localized at 
a Ceo or aromatic hydrocarbon molecule located at 
and Rij=Ri— Rj. The parameters f„„i(Rij) represent 
an onsite energy (R^ =0) and hopping integrals (R^jT^O), 
which are described with the translational symmetry as 

tnm(R) = {4>nR\'HKs\4'mo), (5) 

where |0nRi)=ai„|O) and Hks is the Kohn-Sham Hamil- 
tonian. 

To evaluate effective interaction parameters [/„„(R) 
and Jra„(R), we calculate the screened Coulomb inter- 
action W{r,r') at the low-frequency limit. We first cal- 
culate non-interacting polarization function x with ex- 
cluding polarization processes within the target bands. 
Note that screening by the target electrons is considered 
when we solve the effective models, so that we have to 
avoid double counting of it when we derive the effective 
models. With the resulting x, the W interaction is cal- 
culated as W=(l — vx) V, where i; is the bare Coulomb 
interaction v{r, r')= ■ 

Once screened Coulomb interaction W{r,r') is calcu- 
lated, the matrix elements of W are obtained as 

U,m(R) = ( (pnR^'mO | | 0riR0mO ) 

= I ^drdr'0:^(r)0^(r)M^(r,r')0:^(rV™o(r') (6) 

and 

-^rjm(R) = ( 'Z'nR^rTiO | 1 0^0 ^rjR ) 

= J jdvdv'(p;^(^)cf>,^m{v,v)4,lJc)Mv). (7) 

For comparison with the cRPA results, we calculate in- 
teraction parameters with different levels of screening. 
One is the unscreened one, i.e., the bare Coulomb inter- 
action, and the other is the fully-screened one where we 
calculate x including the target-band screening. To dis- 
tinguish these from cRPA, we denote them as 'bare' and 
'fRPA'. 



III. RESULTS 

A. Calculation conditions 

We performed DFT band calculation with Tokyo 
Ab initio Program Package^ based on the pseudopo- 
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TABLE I: Basic property of fee and A15 alkali-doped Ceo 
eompounds; the lattice parameter a, corresponding Ceo^~ vol- 
ume in solid, and measured superconducting transition tem- 
perature Tc or the Neel temperature Tjv. For fee CssCeo, the 
three samples are specified with the Cg,^ volume and corre- 
sponds to those in the superconducting phase with maximum 
Tc (VgQ^'^), in the vicinity of the metal-insulator transition 
(Vmit), and in the anti-ferromagnetic insulating phase (Vafi), 
respectively. For A15 structure, we also list another sample 
with a higher pressure, for which Tc is lowered than that of 
V^q"'^ , and is abbreviated to V^sc^*^^- 







a 

(A) 


Volume/Ceo^ 

(A3) 


Pressure Tc(rjv) 
(kbar) (K) 


Ref. 




K 


14.240 


722 





19 


49 


fee 
A3 Ceo 


Rb 
Cs(Vmit) 


14.420 
14.500 
14.640 


750 
762 
784 



7 
2 


29 
35 
26 


49 
4, 




Cs(Vafi) 


14.762 


804 





(2.2) 








11.450 


751 


15 


35 




A15 


, ropt.P 


11.570 


774 


7 


38 




CssCeo 


Vmit 


11.650 


791 


3 


32 






Vafi 


11.783 


818 





(46) 





TABLE II: Lattice parameters for pristine solid picene, 
coronene, and phenanthrene and superconducting transition 
temperature Tc observed for doped systems. 





a 

(A) 


b 

(A) 


c 

(A) 


/? 
(°) 


Tc 
(K) 


Ref. 


picene 


8.480 


6.154 


13.515 


90.46 


18,7 




coronene 


16.094 


4.690 


10.049 


110.79 


15 




phenanthrene 


8.453 


6.175 


9.477 


98.28 
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tential plus plane-wave framework. We used the 
generalized-gradient approximation (GGA) exchange- 
correlation functional with the parameterization of 
Pcrdew-Burkc-Ernzcrho^ and the TrouUiar-Martins 
norm-conserving pseudopotentials^ in the Kleinman- 
Bylander representation4i The pseudopotentials for al- 
kali metals, K, Rb, and Cs were supplemented with par- 
tial core corrections^ The cutoff energies for wavefunc- 
tions and charge densities were set to 36 Ry and 144 Ry, 
respectively, and we employed 5x5x5 fc-point sampling. 
We confirmed that this condition ensures well converged 
results. 

The DFT calculations were performed for the follow- 
ing twelve materials: fee KaCgo, fee RbaCgo, fee CsaCgo 
with three different lattice parameters, A15 CssCgo with 
four different lattice parameters, doped solid picene, 
coronene, and phenanthrene. The lattice parameters 
were taken from the experiments and internal coordi- 
nates were optimized.— In fee A3 Ceo, the disorder of the 
orientation of Cgo molecules was ignored, so the crystal 



symmetry is lowered from Fm3m to Fm3. 

Before presenting the computational results, we sum- 
marize the basic properties of the compounds studied 
in the present paper. Table |T] lists experimental val- 
ues for the Cgo eompounds, including the lattice con- 
stant a, the volume per Cgo"^" in solid,— applied pres- 
sure, and measured superconducting transition temper- 
ature Tc or the Neel temperature Tjv- The a value 
and/or Cgo'^^ volume can be controlled by the chemi- 
cal and external pressures. In this table, the samples 
are arranged in the order of the increase of the lat- 
tice constant. For fee A3C60, Tc first increases and 
reaches the maximum (35 K) around a=14.500 A. Then, 
it decreases down to Tc ^ 25 K where the system ex- 
periences the metal-insulator transition (MIT) and be- 
comes an antiferromagnetic insulator (AFI), for which 
the Neel temperature Tn is around 2.2 K. A similar be- 
havior is observed in the A15 system, while T/v is signif- 
icantly higher (46 K). This is because the A15 struc- 
ture is bipartite and therefore less frustrated 1^ Here- 
after, we label the nine Ceo compounds as fcc-K, fcc- 
Rb, fcc-Cs(yg°^*-^), fcc-Cs(yMiT), fcc-Cs(T/AFi), A15- 
Cs(K3c'^''^), A15-Cs(l/s°P' '^), Al5-Cs(yMiT), and A15- 
Cs(1/afi). 

Table [U shows the experimental lattice parameters 
for undoped solid picene, coronene, and phenanthrene 
and Tc observed for doped systems. For doped solid 
picene, two different Tc (18 K or 7 K) have been ob- 
served depending on the preparation conditions.— The 
superconductivity appears when the system is doped, 
but the details of the crystal structures in the super- 
conducting phases have not been determined. Thus in 
the present study, the band calculations for aromatic 
compounds were performed for the artificially charged 
system where three negative charges per one hydrocar- 
bon molecule were doped with a uniform compensat- 
ing positive background charge. Hereafter, doped solid 
picene, coronene, and phenanthrene are referred to as 
solid picene"^", coronene"^", and phenanthrene"^ ~, respec- 
tively. 

B. Band structure and density of states 

Figure [T] shows our calculated GGA band structures 
for the fee A3C60 (upper 5 panels), A15 CssCgo (mid- 
dle 4 panels) , and aromatic compounds (lower 3 panels) . 
These compounds have common features in their band 
structure; i.e., we see narrow bands near the Fermi level 
separated from other bands, being preferable when we 
choose the target bands to construct an effective model. 
In the Cgo compounds, there are threefold degenerated 
states, which form the so-called 'lu band' near the Fermi 
level, and we construct effective models for these bands. 
For aromatic compounds, the target bands are made from 
the lowest two unoccupied molecular orbitals (LUMO 
and LUMO-l-1) in an isolated molecule .i^^i^^i^ It should 
be noted that unoccupied bands lie above the target 
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fcc-K fcc-Rb fcc-Cs(CP'-^) fcc-Cs(K,.J fcc-Cs(F^j) 
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A15-Cs(Fh'gh/') AlS-CsCK^f-^) A15-Cs(F^„) A15-Cs(F^,,) 




^XRMXr RMrXRMXr RMrXRMXr RMrXRMXr RMr 




rx AY rzD Eczrx ay tzd ec ztxaytzd ec z 



FIG. 1: (Color online) Calculated ab intio electronic band structure of fcc-K, fcc-Rb, fcc-Cs( Vg°(S''^), fcc-Cs(l/MiT), fcc-Cs(VAFi), 
A15-Cs(l/s°^'--^), A15-Cs(l/g''c^.«^''^), A15-CsCI/mit), A15-Cs(1/afi), solid picene , solid coronene , and solid phenanthrene . 
In the case of aromatic compounds with monoclinic structure, the horizontal axis is labeled by the special points in the Brillouin 
zone with T, X, A, Y, Z, D, E, and C, respectively, corresponding to (0, 0, 0), (1/2, 0, 0), (1/2, 1/2, 0), (0, 1/2, 0), (0, 0, 1/2), 
(1/2, 0, 1/2), (1/2, 1/2, 1/2) and (1/2, 1/2, 1/2) in units of (a*,b*,c*). The interpolated band dispersion with the derived 
tight binding Hamiltonian is depicted as blue dashed lines. 
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bands more densely in the order of solid picene^" , solid 
coronene'^", and solid phenanthrene^". Since conduction 
bands can generate stronger screening when they reside 
closer to the target bands, we expect a weak repulsive 
interaction in solid picene'^" compared to the other two. 

We show in Fig. [2] the calculated density of states 
(DOS) of the tiu band for fee AaCgo (a) and A15 CssCgo 
(b). For both fee or A15, the bandwidth W monoton- 
ically increases as decreasing the lattice constant, but 
the DOS profile does not change drastically. We list the 
values of W in table HTH The bandwidth of A15 (-0.6 
eV) tends to be larger than that of fee (—0.4 eV), which 
is due to the difference in the inter-Ceo contact, i.e., 
'hexagon-to-hexagon' configuration for A15 and 'bond- 
to-bond' one for fcci^ It was found that the bandwidths 
of the aromatic compounds are nearly 0.5 eV (see ta- 
ble Ell). 
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FIG. 2: (Color online) (a) Our calculated density of states 
(DOS) for lu band of fcc-K (red), fcc-Rb (green), fcc- 
Cs(Vg°Q*'^) (blue), fcc-Cs(VMiT) (purple), and fcc-Cs(VAFi) 
(light blue). (b) DOS for fu band of A15-Cs(l/sc«^''^) 
(red), A15-Cs(ys°<S'-^) (green), A15-Cs(1/mit) (blue), and 
A15-Cs(Vafi) (purple). 




FIG. 3: (Color online) Isosurface of our calculated p^-like 
maximally localized Wannier orbital of A15-Cs(Vafi) viewed 
along the (a) x axis, (b) y axis, (c) z axes, drawn by VESTA.— 
The red surfaces indicate positive isosurface and the blue sur- 
faces indicate negative isosurface. 



C. Maximally localized Wannier orbitals 

Figure [3] shows a contour plot of one of MLWOs for 
the tiu bands of A15-Cs(Vafi). The results of other Ceo 
compounds are almost the same. From this figure we see 
that the resulting Wannier orbital is well localized at the 
single Ceo molecule. In this plot, we displayed the same 
orbital along the three directions; panels (a), (b), and 
(c) correspond to the view along the a;, y, and z axis, 
respectively. We see a node in the center of this orbital 
for (b) and (c), thus this orbital has pj;-like symmetry. 
Note that the view along the y axis is not identical to the 
view along the z axis, which is in contrast with the case of 
atomic p orbitals. We note that the other two Py- and Pz- 
like Wannier orbitals are symmetrically equivalent to the 
presented Pa:-like orbital. We also note that the weight of 
the Wannier orbitals concentrates in the vicinity of the 
cage of the Ceo molecule and there is little weight inside 
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TABLE III: Calculated bandwidth W of the target band and spatial Wannier spread O for twelve materials: fcc-K, fcc-Rb, fcc- 
Cs(Fs°c*'^). fcc-Cs(yMiT), fcc-Cs(VAFi), A15-Cs(K3c'^''^), A15-Cs(Vs°<S*-^), A15-Cs(1^mit), A15-Cs(T/afi), solid picene^', solid 
coronene^", and solid phenanthrene^" . For the aromatic compounds, the two value of Q, axe Usted; the left is the 'lower-level' 
orbital and the right is 'higher-level' one. Units are given in meV for W and A for f2. 









fee AgCeo 






A15 CsaCeo 




aromatic compounds 




K 


Rb 


Cs(K?c"") Cs(14iit) 


Cs(\/afi) 




Vgo"' Vmit 


Vafi 


picene'^^ coronene'''^ phenanthrene''^ 


w 


502 


454 


427 379 


341 


740 


659 614 


535 


477 447 505 


n 


4.28 


4.21 


4.19 4.14 


4.10 


4.27 


4.20 4.16 


4.12 


4.08, 4.13 3.64, 3.67 3.20, 3.08 




IttJI'^ is ^268 A'^. The value is compared with the Cgo'^ 
volume listed in table|I] (^720-820 A"^), clearly indicating 
well-localized nature of MLWO on the single molecule. 
We see that the Wannier spread has a weak positive cor- 
relation with the bandwidth W. In the aromatic com- 
pounds, the molecular size itself is different from each 
other, which result in the appreciable difference in il. 
Note that, O has no clear correlation with W. 



FIG. 4: (Color online) Isosurface of maximally localized Wan- 
nier orbitals of solid phenanthrene^" with (a) lower onsite en- 
ergy and (b) higher onsite energy, drawn by VESTA.— The 
red surfaces indicate positive isosurface and the blue surfaces 
indicate negative isosurface. 



it. 

We next show in Fig.|3]a contour plot of two MLWOs of 
solid phenanthrene^". In the aromatic compounds, the 
two basis orbitals of the effective model are not symmet- 
rically equivalent and therefore we specify these orbitals 
as 'lower' and 'higher' orbitals in terms of the onsite level 
of MLWOs. The lower and higher orbitals are shown in 
the panels (a) and (b), respectively. We again see the re- 
sulting orbitals are well localized at the single molecules. 
The MLWOs of solid picene'^" and coronene'^" are sim- 
ilar to those of undoped systems calculated in Refs. [isl 
andlH 

We list in table IIIll our calculated spatial spread r2„ of 
MLWO for the twelve materials, where f2„ is defined as 

= \/ (^nok^l^no) " | (0nO |r | ^no) | ^ ■ (8) 

In the Cgo compounds, the calculated Wannier spread 
is roughly 4 A and thus the estimated effective volume 



D. Transfer integrals 

Let us move on to transfer integrals. For the Cgo com- 
pounds, the band dispersion of the target band was found 
to be well reproduced only with nearest neighbor (NN) 
and next nearest neighbor (NNN) transfers. The orbital 
index, 1, 2, and 3 denote Px-, Py-, and p^-like orbitals, 
respectively. Onsite energies for three MLWOs are set 
to zero. From now on, 'site' means one molecule and 
the coordinate of site R is defined as the center of the 
molecule. The transfer integrals f„„i(R) are represented 
as 3 X 3 matrix. In fee (A15) structure, there are 12 (8) 
NN sites and 6 (6) NNN sites per site. From transfers 
to the specific site, other transfers to the equivalent sites 
are reproduced by proper symmetry operations. As a 
representative NN site, we choose R=(i?i,, Ry, i?z)=(0.5, 
0.5, 0.0) and (0.5, 0.5, 0.5) for fee and A15 structure, re- 
spectively, where the coordinate is based on conventional 
cell. The transfer matrix to this site is 

FiF2 0\ Ml ^(3) ^(2)\ 

F2F3 and ^(2) Ai ^(3) (9) 
0^4/ \^2(3) ^3(2) A, J 

for the fee and A15 structure, respectively. In A15, the 
two Geo molecules in the unit cell (denoted as A- and 
B-site) are not equivalent in terms of their orientations. 
So, in the matrix (|9]), we show both the transfers from A- 
site to B-site and from B-site to A-site (in parentheses). 
We choose R=(l, 0, 0) for a representative NNN site for 
both structure, then the transfer matrix is written as 

F5 \ [Ai \ 
and ^5(6) (10) 

OF7/ V ^6(5)/ 

for the fee and A15 structure, respectively. The transfer 
matrix for A15 represents the A-A transfer and the B-B 
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TABLE IV: Hopping parameters for fee A3 Ceo in Eqs. ^ 
and (|10p . Units are given in lO^^eV. 





Fi 


Fa 


Fa 


F4 


Fs 


Fe 


F7 


fcc-K 


-40 


-339 


421 


-187 


-94 


-14 


-2 


fcc-Rb 


-16 


-306 


392 


-159 


-75 


-8 


15 


icc-CsiVsT'') 


26 


-299 


372 


-120 


-60 


-3 


36 


fcC-Cs(VMIT) 


15 


-267 


332 


-104 


-40 


1 


30 


fce-CsCl/AFi) 


13 


-241 


302 


-94 


-33 


1 


24 



For A15 structure, we have 



13(2) 



^2(3) 

Ai 

-^3(2) 



A. 



Ai 

A.. 



2(3) 



3(2) 

2(3) "^3(2) 



i3(2) 
2(3) 



-A. 



Ai 
A^ 



Ai 



2(3) ^3(2) 



-A 



-A; 

~A 
Ai 

l2(3) ^3(2) 



3(2) 
2(3) 



2(3) 



A 



Ai 



3(2) 
2(3) 



for R = (0.5,0.5,-0.5) 



for R = (0.5,-0.5,0.5) 



forR= (0.5,-0.5,-0.5). 



TABLE V: Hopping parameters for A15 CsaCeo in Eqs. @ 
and (|10|) . Units are given in 10~'*eV. 



Ai A2 Az Ai 



A. 



Aa 



AT5^ci(V^^f^ 
A15-Cs(K?c*"^) 
A15-Cs(Vmit) 
A15-Cs(1/afi) 



-297 448 67 74 -105 -289 

-262 400 61 74 -97 -239 

-239 371 57 76 -89 -212 

-206 329 53 73 -79 -180 



one (in parentheses) . Table IIVI shows the value of the 
parameters from Fi to Fj for fee Ceo compounds. We 
note that Fq^F-; is due to the lowering the symmetry 
from Fm3m to Fm3. Table IVl shows the value of param- 
eters from Ai to Aq for A15 CsaCgo. For both systems, 
the values of hopping parameters decrease as the lattice 
parameters increases. 

We next describe the procedure for obtaining transfers 
to the other NN sites or NNN sites. First we consider 
the NN case. For fee structure, the transfer matrices to 
other five NN sites arc written with F1-F4 as follows: 



for R = (0.0,0.5,0.5) 














^1 


F2 





F2 


F3 


Fi 











^1 


-F2 





-F2 


F3 


^^3 





F2 





Fi 





F2 







F3 





-F2 





Fi 





F2 







Fi 




o' 


F2 


^^3 











Fi 



forR= (0.0,0.5,-0.5) 



for R = (0.5,0.0,0.5) 



for R (-0.5,0.0,0.5) 



forR= (0.5,-0.5,0.0). 



The remaining transfers to the NN sites are reproduced 
by the relation t„m(R-) = ijim(-R)- 

Similarly, the transfers to the other NNN sites are de- 
scribed as 



for R = (0,1,0) 




forR= (0,0,1) 



for fee (A15) structure and the remaining NNN transfers 
are generated according to t„m(R) = inrn(— R)- 

Using these NN and NNN transfer parameters, we con- 
struct the transfer part Ht in Eq. of the effective 
model. The band dispersion for the Cgo compounds cal- 
culated from the resulting Ht is depicted as blue dashed 
lines in Fig. [1] from which we see that the original GGA 
band dispersion is satisfactorily reproduced. 

For the aromatic compounds, since there is no simple 
symmetry operation, their transfers are difficult to show 
concisely.— For the aromatic compounds, we describe 
only some characteristic features of the transfers. The 
aromatic compounds are regarded as the stacking layered 
systems, so we expect 2-dimensional hopping structure. 
However, in the present transfer analysis, we found that 
the anisotropy of the transfers is not so simple. Table IVll 
compares the maximum absolute value of the intralayer 
transfers (t||) with that of the interlayer transfers {t±), 
where the intralayer is defined as the ab plane. The in- 
tralayer transfers are further decomposed in the three 
directions and compared with each other (see Fig. [5] for 
the definition of the three directions). The anisotropy 
j-^max^^max-j -g ^^^^ appreciable for solid picene"^", es- 
timated as 20/59~0.34, and phenanthrene^", as '-.^0.49. 
In contrast, the anisotropy of coronene'^" is significant 
and is ~0.16. In the case of coronene'^", the intralayer 
anisotropy is even strong tp^'^/f||^^^=7/87~0.08; this sys- 
tem is almost quasi-one-dimensional chain along the b 
axis. We note that the original GGA band dispersion is 
well reproduced by short-range transfer hoppings (Fig. 
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FIG. 5: (Color online) Schematic picture of intralayer transfer 
t|| along the three directions in aromatic compounds. The 
ellipses indicate molecules. 



TABLE VI: Comparison between the maximum absolute val- 
ues of intralayer transfer ty and interlayer transfer t± for aro- 
matic compounds. For the three directions in ty, see Fig. (5] 
Units are given in meV. 





f max 

111 


.max 

IP 


,max 

113 


I max 


solid picene'^" 


48 


39 


59 


20 


solid coronene^" 


7 


87 


7 


14 


solid phenanthrene^" 


49 


32 


73 


36 



E. Effective interaction parameters 

We performed RPA calculations to evaluate the 
screened Coulomb interaction W{r, r') in Eqs. ([6]) and 
([7|), where the dielectric function was expanded in plane 
waves with an energy cutoff 7.5 Ry for fee A3C60 and 
aromatic compounds and 5.0 Ry for A15 CsaCgo- The to- 
tal number of bands considered in the polarization func- 
tion was set to 335 (120 occupied-|-3 target-l-212 unoc- 
cupied) for fee A3C60, 670 (240 occupied-l-6 target-|-424 
unoccupied) for A15 CsaCeo, 310 (102 occupied-|-4 tar- 
get -1-204 unoccupied) for sohd picene'^", 315 (108 oc- 
cupied-l-4 target -1-203 unoccupied) for solid coronene'^", 
and 270 (66 occupied-|-4 target -1-200 unoccupied) for solid 
phenanthrene'^" . The Brillouin-zone integral on wavevec- 
tor was evaluated by generalized tetrahedron method/^ 
A problem due to the singularity of long-wavelength- 
limit Coulomb interaction in the evaluation of the Wan- 
nier matrix elements, t/„m(R) in Eq. ([5]) and J„m(R)Jji 
Eq. ([7]), was treated in the manner described in Ref. ISOl . 

The onsite interactions are specified by U =C/„„(0) and 
U'~Unm{0) and J=J„,„(0) for n^m. In the case of Cgo 
compounds, J7, fJ', and J take only one value according 
to the symmetry. For aromatic compounds, U is different 
for two orbitals, so we present two values. We also denote 
the Coulomb repulsion between the neighboring sites as 
V. 



Table IVIII shows our calculated interaction parame- 
ters [/', J, and V with three screening levels ('bare', 
'cRPA', and 'fRPA'). We see that the value of the 
Coulomb repulsion decreases as the screening processes 
increases. In the Cgo compounds, the bare value is '^3.4 
eV and after considering the screening by cRPA, the 
value is reduced to '^l eV. By taking account of the 
intra-target-band screening by fRPA, the value is further 
reduced to ~0.1 eV. It should be noted here that the ma- 
terial dependence of the bare values in fee or A15 is small; 
for example, 3.27 eV for fcc-K and 3.37 for fcc-Cs(VAFi)- 
The difference is nearly 3 %. This difference is ascribed 
to the difference in the spatial spread of MLWOs (see 
Table IIII[) . On the other hand, the material difference 
in the cRPA values is beyond 20 %; 0.82 eV for fcc-K 
and 1.07 for fcc-Cs(VAFi)- Indeed, this appreciable dif- 
ference originates from the difference in the macroscopic 
dielectric constant defined as 



^cRPA 



lim lim 

Q-i-O w->0 ,cRPA 



1 



(q,w) 



(11) 



with a; being frequency and Q=q-f-G, where q is a wave 
vector in the first Brillouin zone and G is a reciprocal 
lattice vector. We list the value in the bottom of the 
table. We see that the material dependence of (^^f^ 
is appreciable as 5.6 for fcc-K and 4.4 for fcc-Cs(VAFi), 
clearly indicating the importance of the screening effect 
in addition to the spatial Wannier spread. 

For the aromatic compounds, differences in both the 
bare interaction and the screening effect contribute to 
the material dependence of the cRPA values; the bare 



interactions are V^^^^ 



r rcoroncnc 
'^barc 



bare 

and after consideration of the cRPA screening, we ob- 

tain C/™"< t^cTp™"- t^S;r""™""- Especially, 
in picene'^ , the dielectric constant is markedly high as 
^--^12,— making the cRPA U value small appreciably. 

We finally remark some points: As for the Ceo com- 
pounds, the equality U'^U~2J holds among effective pa- 
rameters. This relationship also roughly holds for the 
aromatic compounds. The present U value of cRPA for 
the Cgo compounds are small compared to the previous 
estimates of U (^l-1.5eV)<^"— For all materials, the 
cRPA Coulomb interaction decays as l/(e^^^^r) with 
r being the distance between the centers of MLWOs, 
while the fRPA Coulomb interaction is limited to be short 
ranged due to the metallic screening (see table IVlI)) . We 
note that the fRPA U gives an opposite trend to the bare 
and cRPA U; for example, in the fee Ceo compounds, the 
fRPA value slightly decreases as proceeding from fcc-K 
to fcc-Cs(VAFi)- This is due to the fact that the Coulomb 
interaction is efficiently screened due to the increase in 
the density of states accompanied by the decrease of 
bandwidth. Wc also found that, in these systems, the 
exchange interaction J are also efficiently screened; i.e., 
'/cRPA/>/baro~0.3. This makcs a clear contrast to the case 
of the inorganic materials as JcRPA/>/bare~0.8 such as 3d 
transition metals,— its oxides SrVOs,— and iron-based 
superconductorsi^2ii£ 
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TABLE VII: U, U' , J, and V with three different screening levels [unscreened (bare), constrained RPA (cRPA), and fully- 
screened RPA (£RPA)] for the twelve compounds: fcc-K, fcc-Rb, icc-CsiVs^"-^), fcc-Cs(\/MiT), fcc-Cs(VAFi), A15-Cs(\/gc'^''^), 
A15-Cs(Vg°Q*'^), A15-Cs(Vmit), A15-Cs(Vafi), solid picene^^, solid coronene''^, and solid phenanthrene''^. For the aromatic 
compounds, the two value of U are presented; the left is the 'lower-level' orbital and the right is 'higher-level' one. For 'bare' 
and 'cRPA' U, U' and V values, the unit is given in eV and J is given by meV. For 'fRPA', the unit is given in meV. In the 
bottom, we present our calculated cRPA macroscopic dielectric constant e^^^ in El- (HU- 



fee AsCgq A15 CS3C60 aromatic eompounds 





K 


Rb 




Cs(Vmit) 


Cs(Vafi) 






ViviIT 


Vafi 


piccnc'^ 


coroncnc'^ 


phcnanthrene"* 




3.27 


3.31 


3.32 


3.35 


3.37 


3.36 


3.39 


3.40 


3.42 


4.43,4.41 


4.64,4.59 


5.05,5.17 


ULr. 


3.08 


3.11 


3.12 


3.15 


3.17 


3.16 


3.18 


3.20 


3.22 


3.55 


4.33 


4.55 




96 


99 


100 


101 


102 


97 


99 


100 


101 


166 


129 


275 


Vbaro 


1.31-1.37 


1.30-1.35 


1.29-1.34 


1.28-1.33 


1.27-1.32 


1.37-1.38 


1.36-1.37 


1.35-1.36 


1.34-1.34 


2.08-2.32 


2.79-2.84 


2.29-2.43 


f-'cRPA 


0.82 


0.92 


0.94 


1.02 


1.07 


0.93 


1.02 


1.07 


1.14 


0.73,0.74 


1.29,1.26 


1.33,1.37 




0.76 


0.85 


0.87 


0.94 


1.00 


0.87 


0.95 


0.99 


1.06 


0.58 


1.15 


1.17 


^^cRPA 


31 


34 


35 


35 


36 


30 


36 


36 


37 


53 


58 


101 


VcRPA 


0.24-0.25 


0.26-0.27 


0.27-0.28 


0.28-0.29 


0.30 


0.30 


0.31 


0.32 


0.34 


0.26 


0.59-0.60 


0.47-0.48 


UmrA 


93 


91 


91 


86 


83 


107 


102 


99 


93 


155,151 


149,120 


166,172 


^tRPA 


41 


39 


39 


35 


32 


50 


45 


42 


37 


51 


53 


60 


-^fRPA 


25 


26 


26 


26 


25 


28 


28 


28 


28 


38 


39 


57 


^4rpa 


1-3 


1-3 


1-3 


1-3 


1-3 


2-3 


2 


2 


1-2 


1-4 


1-4 


2 


^cHPA 


5.6 


5.1 


4.9 


4.6 


4.4 


4.7 


4.4 


4.3 


4.1 


12.0 


5.5 


6.3 



IV. DISCUSSIONS 
A. Material dependence of effective parameters 

Let us move on to the comparison of the effective inter- 
action parameters among the twelve compounds. Figure 
[6]summarizes the results of the cRPA calculation; the on- 
site Coulomb repulsion [7 averaged over MLWOs derived 
from the target band, the onsite exchange interaction J, 
the offsite interaction V averaged over nearest-neighbor 
sites, and the ratio {U — V)/W which measures the cor- 
relation strength in the system. Note that the net inter- 
action is estimated as U — V, based on the analysis of the 
extended Hubbard model. 

As for the Cgo systems, we see that U has apprecia- 
ble material dependence, ranging from 0.8 eV to 1.1 eV 
[Fig. mja)]. This is ascribed to the differences in the 
size of Wannier orbitals and dielectric screening (see Sec. 
IIIIEp . On the other hand, the material dependence of J 
is weak and the value itself is negligibly small as ~0.03 
eV [Fig. Elb)]. In general, small J favors low-spin states, 
as is observed in experiments.—"—'^ It is interesting to 
note that there is a proposal that the Jahn- Teller cou- 
pling dominates over the Hund's rule coupling J, and 
induces superconductivity with the help of sufficiently 
large Ui^ Compared to J, we found that V is substan- 
tially large as ~0.3 eV, being as large as ~25% of [Fig. 
[S]Jc)]. As for W, which measures kinetic energy, we ob- 
served a decreasing trend as the lattice constant increases 
[Fig. m^d)]. We also note that the energy scale for A15 
CsaCgo is larger than that of fee A3 Ceo as mentioned in 
Sec. lIIIB] The derived correlation strength of {JJ — V) /W 
exhibits a rather simple monotonic increasing behavior 
[Fig. ini[e)] , with the lattice constant increase. The pre- 
sented [U — V)/W^\ indicates that Ceo superconductors 
are categorized as strongly correlated electron systems. 

For the aromatic superconductors, we found that the 
energy scale of U is similar to that of the Ceo supercon- 



ductors [Fig. [nja)]. On the other hand, it is interesting 
to note that the aromatic superconductors tend to have 
larger J and V [Figs. |6l[b) and (c)]. We can also see that 
the material dependence of the interaction parameters 
among the aromatic superconductors is also more signif- 
icant, since the size and shape of the aromatic molecules 
are quite different from each other. As for W ^ they are 
similar for the aromatic and Ceo superconductors [Fig. 
[n{d)]. We found that aromatic superconductors are also 
in strongly correlated regime as Ceo ones, based on the 
analysis of the correlation strength [Fig. |6{e)]. 



B. Relation between electronic correlation and 
superconductivity 

Next, let us discuss the relation between electronic cor- 
relation and superconductivity for the Cgo and aromatic 
superconductors. In Figs. [7] (a) and (b), we plot the su- 
perconducting transition temperature Tc and the Ncel 
temperature T/v as a function of volume occupied per 
fuUeride anion (Vb) for fee AaCgo and A15 CssCeo, re- 
spectively (see also table |T| . To see the relation between 
the electron correlation and the superconductivity we su- 
perpose a plot of (C7 — V)/W on the phase diagram. We 
see that while {U — V) /W and Tc have a positive correla- 
tion up to Vb'^760-770 A'^, for larger Vq, electron correla- 
tion becomes fatal for superconductivity and the system 
eventually becomes an insulator. We note that the criti- 
cal value of {U — V)/W for the MIT sample is larger for 
fee A3C60 (~1.9) than A15 CsaCeo (^1-2). As discussed 
in Ref. l56l . it is important to consider the influence of 
lattice and orbital structure on MITjiii^ 

In Fig. [TJc), we plot and {U - V)/W for the 
three aromatic superconductors, which shows a nega- 
tive correlation. Therefore, it seems that electronic 
correlation does not favor superconductivity in these 
aromatic superconductors. Recently, doped 1,2:8,9- 
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fee A3C,„ 



A15 Cs,Q„ 

3 60 



aromatic 
compounds 



FIG. 6: (Color online) Material dependence of (a) the average of the onsite effective Coulomb repulsion U , (b) the onsite 
effective exchange interaction J, (c) the average of the offsite effective Coulomb repulsion between neighboring sites V , (d) the 
bandwidth of the target band W, and (e) the correlation strength [U — V)/W. 
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(a) fccA3C6o (b)A15 CsaQp (c) aromatic compounds 

— ^ ' ' ' • ■ — ■ — ■ — • — 1 p — ' — ■ — ' — • — ' — ■ — ' — T=rn2 <) I I ^-^ 

i r I hl^ l hhll ¥ - IP P 

720 740 760 780 \800 740 760 780 800 820 " • 3. 1 i. 1 

„ 3 , / X 3^ \ A T-T ^ . , ,1^^ picene' coronene^- phenanthrene^- 

Ceo volume (A^) AFI Cgo^" volume (A^) 

FIG. 7: (Color online) Relation between the experimental curve of superconducting or magnetic transition temperature (Tc, T]v) 
as a function of the Cg^ volume and the estimated correlation strength {U — V)/W (vertical bar): (a) fee A3C60 and (b) A15 
CsaCeo- For aromatic compounds (c), the measured Tc in table |II] and the calculated correlation strength are compared, 
where picene^~ = (C22Hi4)^~, coronene^~ = (C24Hi2)^~, and phenanthrene"^~ = (Ci4Hio)"'~. For the panels (a) and (b), the 
experimental phase diagram were taken from Ref. Ufor fee and Ref. for A15. 



dibenzopentacene was found to have a quite high 
K. Since l,2:8,9-dibenzopentacene is a bigger molecule 
than picene, coronene, and phenanthrene, the former 
interaction is expected to be small compared to the 
latter ones, reflecting the large Wannier spread of the 
l,2:8,9-dibenzopentacene molecule. If there is no dras- 
tic change in the bandwidth W, which is probable in 
terms of the tendency shown in Fig. [BJd), the weakest 
electronic correlation will be realized in doped 1,2:8,9- 
dibenzopentacene. This trend is consistent with Fig. 

mc). 

Regarding the role of electronic correlation in the Ceo 
and aromatic superconductors, there are two possibili- 
ties: The pairing mechanism in these compounds has a 
common root or these superconductors have completely 
different pairing glues. If we assume that the aromatic su- 
perconductors reside in the vicinity of the border between 
the superconducting and insulating phases, the first sce- 
nario is (at least partially) explicable to the behavior in 
Fig. [T] On the other hand, in the second scenario, the 
electronic correlation enhances superconductivity for the 
Ceo compounds and inversely suppresses that for the aro- 
matic compounds. In order to clarify this issue, experi- 
mental studies to determine the phase diagram for aro- 
matic superconductors against temperature and volume 
occupied per anion are highly desiredi^ Theoretically, 
microscopic calculations considering both electronic cor- 
relation and electron-lattice coupling are needed, which 
will be an interesting future problem. 

V. SUMMARY 

To provide insight into the role of electronic correla- 
tion in Ceo and aromatic superconductors, we derived 
effective models for wide range of the examples; fcc- 
K, fcc-Rb, fcc-Cs(yg°P*-^), fcc-Cs(yMiT), fcc-Cs(VAFi), 
A15-Cs(K5cS^^), Al5-Cs(Vg°^* -^), A15-Cs(Vafi), solid 



picene'^", coronene^^, and phenanthrene'^^. To define 
the basis orbital of the effective model, we constructed 
MLWOs of isolated bands around the Fermi level. Trans- 
fer parameters are derived by evaluating the matrix ele- 
ments of the Kohn-Sham Hamiltonian between the ML- 
WOs. The low-energy electronic structures of the Ceo 
compounds are highly symmetric and isotropic, so that 
the original GGA band is reproduced with only 6 or 7 pa- 
rameters. On the other hand, the aromatic compounds 
have quite anisotropic electronic structure. 

To quantify the strength of electronic correlation in 
these compounds, we estimated the effective interaction 
parameters such as U, V and J, by means of the cRPA 
method. It was found that, in addition to the appreciable 
reduction of the diagonal part of the Coulomb interac- 
tion (U and V), the off-diagonal part J is also efficiently 
screened. Interestingly, all the Ceo and aromatic super- 
conductors studied in the present work have a similar en- 
ergy scale for the bandwidth and interaction parameters: 
eV, f/-l eV, J-0.05 eV, F~0.3 eV. This param- 
eter range suggests that these compounds arc a strongly 
correlated electron system. However, after examination 
of the material dependence, we found that a clear dif- 
ference between the Ceo and aromatic compounds in the 
relation between electronic correlation strength and T^, 
i.e., a positive correlation in the Ceo system and a nega- 
tive correlation in the aromatic system. 

In the present study, we focused on the derivation for 
the electronic part of the effective model. For deep under- 
standing of the low-energy physics for the carbon-based 
materials, however, the derivation of the electron-phonon 
interaction part is also imperative. The derivation for 
this part includes subtle problems on the definition of 
the basis for the phonon mode (Refs. |63 andlosl) and/or 
the exclusion of the double counting of the screening of 
the low-energy degree of freedoms, which needs the fu- 
ture studies. 
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